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Abstract 



O The construction of the master T-operator recently suggested in [1] is applied to 

integrable vertex models and associated quantum spin chains with trigonometric R- 
matrices. The master T-operator is a generating function for commuting transfer 
matrices of integrable vertex models depending on infinitely many parameters. 
At the same time it turns out to be the tau-function of an integrable hierarchy 
of classical soliton equations in the sense that it satisfies the the same bilinear 
Hirota equations. The class of solutions of the Hirota equations that correspond 
CNj ■ to eigenvalues of the master T-operator is characterized and its relation to the 

^""^ , classical Ruijsenaars-Schneider system of particles is discussed. 



1 Introduction 

(N 



The master T-operator was recently introduced in [T]. It is a generating function for 
commuting transfer matrices of integrable vertex models and associated quantum spin 



chains which unifies the transfer matrices on all levels of the nested Bethe ansatz and 
Baxter's Q-operators in one commuting family. It was also proven in [I] that the master 
T-operator, as a function of infinitely many auxiliary parameters (one of which being the 
usual spectral parameter), satisfies the same hierarchy of bilinear Hirota equations as the 
classical r-function does. Since the operator- valued generating functions commute for all 
values of the auxiliary parameters, there is no problem with their ordering in the bilinear 
equations. 

A similarity between quantum transfer matrices and classical r-functions was first 
pointed out in [5] (see also [3]), where a discrete integrable dynamics in the space of 
commuting integrals of motion of a quantum integrable model was introduced. This 
classical dynamics was identified with the discrete 3-term Hirota equation with special 
boundary conditions. The diagonalization of transfer matrices by means of the nested 
Bethe ansatz technique was shown to be equivalent to an "undressing" chain of Backlund 
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transformations for the discrete Hirota equation. Later this approach was extended 
to supersymmetric integrable models |4j. An essential further step was made in the 
important paper [5], where an operator realization of the Backlund flow describing the 
"undressing" process was constructed for generalized quantum spin chains with rational 
GL(iV)-invariant R- matrices. In fact the master T-operator was already used implicitly 
in that construction. A more explicit and more general definition was given in [I]. 

In this paper we review the construction of [1J trying to avoid technical details. Here 
we deal with the class of integrable lattice vertex models of statistical mechanics with 
trigonometric i?-matrices. The main claim is that the master T-operator for these models 
is a r-function of the classical MKP hierarchy. 

We also characterize the class of solutions of the Hirota equations that correspond 
to eigenvalues of the master T-operator and make explicit the close connection with the 
classical Ruijsenaars-Schneider system of particles [6] which emerges as the dynamical 
system for zeros of the (eigenvalues of) the master T-operator. In an equivalent way, 
the connection emerges from the Baker- Akhiezer function for the Ruijsenaars-Schneider 
system which generates the algebra of commuting operators (transfer matrices) for the 
vertex model (the Bethe algebra). As is well known, the Ruijsenaars-Schneider model 
contains the Calogero-Moser system of particles as a limiting case. In this connection 
let us note that a similar relation between the quantum Gaudin model (which can be 
regarded as a degeneration of quantum spin chains or vertex models with rational R- 
matrices) and classical Calogero-Moser system was found in [7] from a different reasoning. 

2 The transfer matrices 

We consider generalized quantum integrable vertex models with trigonometric i?-matrix. 
The simplest i?-matrix is the operator in (g> C N of the form 

N 

R(u) = ( e T(«+D _ e -7(*+i)) £ 6aa g, ebb + (e 7« _ e -7«) J2 e aa ® e bb 

0=1 Ka^b<N 

(2.1) 

+ (e 7_ e -7) J- e^ b - a ^ u e ab ®e ba . 

l<a^b<N 

Here u G C is the spectral parameter and e ab denotes the N x N matrix with 1 in 
position (a, b) and elesewhere. The deformation (anisotropy) parameter 7 is assumed 
to be such that q = e 7 is not a root of unity. Following the tradition, we call this R- 
matrix trigonometric although the coefficients are hyperbolic functions of u like smh^u. 
Let Vi = C N be copies of the space C , then by R%j(u) denote the i?-matrix acting in 
Vi <8> Vj. The i?-matrix (12.11) satisfies the Yang-Baxter equation 

Ri 2 (tti - u 2 )Rn(u 1 - u 3 )R 23 (m2 - u 3 ) = R 23 (m 2 - u 3 )Ri 3 (ui - u 3 )R 12 (u 1 - u 2 ), (2.2) 

where the both sides are operators in V\ ®V 2 ®V 3 . For any diagonal N x N matrix g 
set gi — g <S> I, g 2 — I ® 9, then the R- matrix commutes with g\g 2 . 

R12H 9i92 = 9\92 Ri 2 (w)- (2.3) 
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This property will be referred to as g-invariance of the i?-matrix. 

Fix a diagonal matrix g = diag(pi,p2> • • • >Piv)- We call it the twist matrix. Below 
we assume that all pi G C are in general position, i.e., Pi/pj ^ e 2in for any i ^ j and any 
integer n. The transfer matrix for the vertex model with twisted boundary conditions 
and with inhomogeneity parameters U{ at each site is defined as 

T(u) = tr (R 10 (u - u 1 )R 20 (u - u 2 ) . ■ .R L o(u - u L ) g^. (2.4) 

The i?-matrices and g are mulitiplied as matrices in the common space Vq (the auxiliary 
space). Trace tro is taken in the auxiliary space. The result is an operator acting in the 
tensor product of vector representations V, = <$f =1 Vj = (C^) 81, (the quantum space). 
Formally, our setting includes also models with higher representations at the sites because 
they can be obtained by "fusing" several vector representations with properly chosen 
parameters Mj. By construction, the operator (12. 4p is a Laurent polynomial in e lu . 

It follows from the Yang-Baxter equation and from the g-invariance of the i?-matrix 
that the transfer matrices for models with the same 7 and g commute for all u and can 
be diagonalized simultaneously. Their diagonalization is the basic problem of the theory 
of vertex models. The standard method is the nested Bethe ansatz technique. 

The full commutative family of operators in the quantum space is in general larger 
than the one generated by coefficients of T(u). The algebraic construction of higher 
commuting transfer matrices essentially relies on representation theory of the g-deformed 
algebras U q (gl(N)) and U q (gl(N)) (see, e.g., [H [EH H21 M ) • 

The algebra U q (gl(N)) has generators L^ b with 1 < a < b < N and L~ b with 1 < b < 
a < N such that L+ a L~ a = L~ a L+ a = 1. Combining them into matrices 

L + = e *b ® L^ b , L - = J2 e *b ® Kb 

a<b a>b 

with U q (gl(N))-vahiedi matrix elements, one can represent the defining relations of the 
algebra in the form 



Ri2L^L^ — L^L^Ri2, Ri2L^L 2 — L 2 L^"Ri2 
with the M-independent i?-matrix 

JV 

R12 = lim (e~ lu R(u)) = ?Ve (M 8e it + V e aa <g> e bb + (q-q^ 1 ) V e ab ® e ba . 

g i" — ^OO 

a=l l<a=£b<N l<a<b<N 

The diagonal elements L^ a can be understood as exponents of the commuting Cartan 
generators h a : 

r± _ „±h a 

^aa — 1 

Let 7t\ be the irreducible finite-dimensional representation of U q (gl(N)) with the 
highest weight A = (Ai, A 2 , • • • , Ajv) such that Aj G Z + , Ai > A 2 > . . . > Ajv. The highest 
weight vector v obeys 

L~ b y = 0, a > b, L" v = q- Xa v. 
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The representation space is generated by repeated action of the generators L\ on 
the highest weight vector and subsequent factorizing (see [TT] for details). These repre- 
sentations are g-deformations of the highest weight finite-dimensional representations of 
U(gl(N)). The highest weights are naturally identified with Young diagrams (partitions) 
A. 

The representation irtu corresponding to the one-box diagram is the vector repre- 
sentation in C . On the Cartan generators h a introduced above it looks exactly like 
for the usual non-deformed algebra gl(N): 7T(i)(h ) = e aa . Given the diagonal matrix 
g = diag(pi,p 2 , • • • ,Pn), we set 

g = Pi 1 P2 2 ---PN N , (2-5) 

then g = 7r (1) (g). 

The iZ-matrix R A (w) acting in C N <g> V w is 

R x (u) = e' u J2 e ab ® tt a (L+ ) - £ e ab ® tt a (L; 6 ). (2.6) 

a<b a>b 

The i?-matrices R A (w), R m (m) are intertwined by a more general i?-matrix R Xfl (u) which 
acts in V {x) ® 

R A ^(Ml - U 2 )R X 3 (U 1 - U 3 )R% 3 (lL 2 - Us) = R 23 (u 2 - W 3 )R A 3 (wi - M 3 )Ri2 («i - u 2 ) (2.7) 

This Yang-Baxter relation generalizes (12.21) . Here space 1 is V^ x \ space 2 is and 
space 3 is C N . The explicit form of R XfM (u) is much more complicated than (12.61) . It can 
be obtained from the universal R- matrix for the quantum affine algebra U q (gl(N)) [15] by 
specifying it to finite-dimensional evaluation representations or by the fusion procedure 
[16l[T71[T3] applied to the fundamental i?-matrix R(u). The g-invariance (12.31) is extended 
to R x »(u) as follows: 

Rl2» 7T A (g)l 7T„(g) 2 = 7T A (g)l 7T M (g) 2 R^( U ). (2.8) 

The higher transfer-matrices, or T-operators, are constructed in a similar way to (12.41) 
by taking trace of the product of R- matrices R x (u — ui) in the auxiliary space 

T x {u) = tr vW (R A (w - nx)R A (n - u 2 ) . . . R x L0 {u - u L ) 7T A (g)). (2.9) 

Here the space with index is the auxiliary space V^ x \ These T-operators act in the 
same quantum space H = C® L . If A = (1) is the 1-box diagram, then definition (12.91) 
coincides with (12.71) . By analogy, we will call g of the form (12. 5p the twist element. 
The Yang-Baxter equation (12.71) and the g-invariance (I2.8P imply that the T-operators 
with the same g commute for all u and A: [T x (u), T^{v)] = 0, and can be diagonalized 
simultaneously. 

An important property of the T-operators defined by (12.91) is that they vanish iden- 
tically if the first column of A is longer than N. 

Set 

L 

H a = Y, e aL e£ = I®_. . <g> I® e aa ® I ® . . . ® I , 
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then the g-invariance implies that [T x (u), H a ] = 0. Therefore, the eigenstates of the 
transfer matrices can be classified according to eigenvalues of the operators H a . Let 

U= U{{M a }) 

Mi,...,M N 

be the decomposition of the quantum space H into the direct sum of eigenspaces for the 
operators H a with the eigenvalues M a G Z+, a = 1, 2, . . . , N, then eigenstates of T x (u) 
lie in the spaces 7i({M a }). Since J2a e aa = H is the unit matrix, J2 a H a — LI® L and thus 

JV 

£M i = L. (2.10) 

i=i 

For the trivial representation (corresponding to the empty Young diagram 0) tt%{L^) = 
1 if a = b and otherwise, n^g) = 1. Formula (12.61) yields _R (m) = e 7 " — e~ 7 ", where 
multiplication by the unity matrix I is implied. Therefore, we can define the T-operator 
for the trivial representation as follows: 

L 

T\u) = 2 L Y[smh^(u- Ui )). (2.11) 

For the one-dimensional representation 7T(iw\ (corresponding to the Young diagram (1^) 
with one column of height N) 7r ( - 1 jv- ) (g ha ) = q for diagonal generators and otherwise, 
7T0(g) = detg. Formula (I2.6p yields R( lN \u) = e l( ~ u+1 ^ — e~ l( - u+1 \ where multiplication 
by the unity matrix I is implied. Therefore, the T-operator for the representation with 
the highest weight (1^) (the quantum determinant of the quantum monodromy matrix) 
is given by: 

L 

T {lN \u) = 2 L det#n sinh (7( M + 1 ~ u i)) = detgT®(u + l). (2.12) 
For general A the T-operator is the Laurent polynomial in e 7n of the similar form: 

L/2 

T\u)= G x k e 2k ^ u . (2.13) 

k=-L/2 

The coefficients G x of the T-operators with fixed g generate the full family of commuting 
operators (the Bethe algebra of the vertex model). 

The operators T x (u) appear to be functionally dependent. They are known to obey 
some functional relations which are given by the Cherednik-Bazhanov-Reshetikhin (CBR) 
determinant formulas [H [9]. These formulas express T x (u) for arbitrary A through the 
transfer matrices T s {u) := T^ s \u) corresponding to 1-row diagrams of length s or through 
the transfer matrices T a {u) := T^iu) corresponding to 1-column diagrams of height a: 

T\u) = ( J] T\u-k)) det T Ai _ l+J (n-j + l), (2.14) 

Ai-1 

T\u) = ( TT T\u + k)) det T x '*- l+j {u+j-l). (2.15) 

Hereafter A' denotes the transposed diagram (with respect to the main diagonal), so that 
A' x is the height of the first column, and is the empty diagram. One can show that 
formulas ( I2.15P follow from (I2.14p and vice versa. 
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3 The master T-operator 



Let t = {ti,t2,ts, . . .} be an infinite set of parameters which we call times because they 
will have the meaning of hierarchical times in the MKP hierarchy. The Schur polynomials 
s\(t) labeled by Young diagrams A can be defined by the determinant formula 

s A (t) = . det ^ f h Xi - i+j (t), (3.1) 

where the polynomials hj are defined with the help of the generating series 

exp(j2 hz k ) = 1 + h x (t)z + h 2 (t)z 2 + ... 

k>l 

It is convenient to put ho(t) = 1, h n (t) = for n < and S0(t) = 1. The functions hj 
are elementary Schur polynomials in the sense that for 1-row diagrams A = (j) with j 
boxes sy)(t) = hj(t). Equivalently, one can define 

s A (t)= det ey_ l+J (t), (3.2) 

i,j=l,...,Ai % 

where the polynomials Cj are defined with the help of the generating series 
exp^-l)*" V fc ) = 1 + ei {t)z + e 2 (t)z 2 + ... 

k>l 

For 1-column diagrams A = (P) with j boxes S(u)(t) = e 3 -(t). Equations ( 13.11) . f ]3 .2f) are 
known as Jacobi-Trudi formulas. It can be proved [18] that the Schur polynomials form 
a basis in the space of symmetric functions of the variables Xi defined by ktf. = J2i x i- 

We note the Cauchy-Littlewood identity 

5> A (t)s A (t') = exp(Y, kt kt' k ), (3-3) 

A k>l 

where the sum is over all Young diagrams including the empty one. Writing it in the 
form 

^s A (y)s A (<9) =exp(Y / ykd t \ 



X fc>l 

la i. 



where d = {d tl , |9 t2 , ^d t3 , . . . } and applying to s M (t), we get: 



s\(8)sJt) 



t=0 



Sxu- (3.4) 



Following PP, we introduce a generating function of the T-operators (the master T- 
operator) depending on the infinite number of parameters t = {t\, t 2 , ■ ■ ■}'■ 

T(u, t) — s\(t)T x (u). (3.5) 

A 

These operators commute for different values of the parameters: \T(u,t), T(u',t')] = 0. 
Since T x (u) = if A^ > N, the sum in (13.51) is actually restricted to diagrams with 
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X[ < N. The T-operators T x (u) can be restored from the master T-operator according 
to the formula 

T\u) = s x {d)T{u,t) 



t=o 



(3.6) 



which follows from (13. 4p . In particular, 

T\u) = T(u,0), (3.7) 

T(u)=T^(u)=d tl Tfat) . (3.8) 

t=o 

Below we use the standard notation 

t ± [z- 1 ] = {ij ± z~\ t 2 ± i f 3 ± i z-\ ...}, (3.9) 

oo 

£(t,z) = ]>>z fc . (3.10) 

k=0 

Eq. (13.61) implies that T(u, 0±[z~ 1 }) is the generating series for T-operators corresponding 
to the 1-row and 1-column diagrams respectively: 

T(u, [z- 1 ]) = £ z- s T s (u), T(u, -fa 1 ]) = £(-*)-"T>)- (3.11) 

s>0 a>0 

As it was proven in pQ, the CBR formulas ( 12. 14ft imply that the master T-operator 
obeys the bilinear identity 

£ e at-f, z) z u-u' T ^ t _ ^-ijj t(m / ; t , + ^-ij) dz = Q (3>12) 

for all u, u', t,t'. The contour C encircles the cut between and oo. By standard ma- 
nipulations [HI [20] , one can derive from (I3.12p the infinite KP and MKP hierarchies of 
differential (in £&'s) and differential- difference (in t^s and u) equations. The variable u 
is the so-called "zero time"; it is naturally included in the extended sequence of times 
to = u, ti, t 2 , . . .. Choosing u, u', t, t' in a special way, one can also derive from (I3.12p the 
following bilinear equations: 



(z 2 - z 3 )T (u, t + fa 1 ]) T (u,t+ fa 1 ] + fa 1 ) 
+ (z 3 - z x )T (u,t+ fa 1 ]) T (u, t + fa 1 } + fa 1 ]) (3.13) 
+ ( Zl - z 2 )T (u, t + fa 1 }) T (u, t + fa 1 } + fa 1 ]) = 0, 



z 2 T (u + 1, t + fa 1 ]) T (u, t + fa 1 }) - Zl T (u + 1, t + fa 1 ]) T (u, t + fa 1 ' 
+ fa - z 2 )T (u + 1, t + fa 1 ) + fa 1 ]) T(u, t) = 0. 



(3.14) 



They are known as Hirota or Hirota-Miwa equations for the r- function [21] [22] . In this 
sense the master T-operator (any of its eigenvalues) is the r-function of the classical 
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MKP hierarchy (see, e.g., [23J). Equation (I3.5P can be regarded as the Schur function 
expansion of the r-function (see also [23]). 

Note that the transformation 



T(u,t) ->■ C(u) exp(j2c k t k ^T(u,t) 



with arbitrary function C(u) and arbitrary constant coefficients c k preserves the space of 
r-functions. Two r-functions are regarded as essentially different if they are not obtained 
from each other by such transformation. 

As a function of u, the master T-operator has the structure similar to ( 12. 13(1 : 



L/2 

T(u,t)= J2 G*(t)e 2fc7W . (3.15) 

k=-L/2 

In particular, the highest and the lowest coefficients G±L/2(t) are easy to calculate. For 
example, the highest coefficient G^^it) is 

L 

G x L/2 (t) = exp(- 7 £ Un) £ tT vw ( L tbi L U 2 • • • L a L b L ) e^ejia • • • e S> 

n=l o.i,bi 



L 



(one should take the first terms from each i?-matrix (12.6(1 ) . Since all matrices L\ h here 
are upper triangular, the trace is equal to that of the product of diagonal matrices with 
the same diagonal elements: 

G£ /2 (t) = exp(- 7 £ u n ) 5> vw (q^+~+^ . ..ft) e« 6i . . . e ™ L 



n=l a, 
L 



£ exp(- 7 E^)tr y( ,)((e^V) hl ...(e^)^) £ e£> ai . . . . 



v aiai ' ' ' 1 

Mi,...,A/ L ' n=l ' 'oi:{Mj} 



where the last sum goes over all sequences of indices a\, . . . such that the number of 
indices equal to j is Mj. It is easy to see that 

V = U p.. 



1 " 



where Pm 1: ...,Mn * s the projector to the subspace %({Mj}). Set y k = — £e 7Mafc p^, then 

^ o=l 

tr^cx, ((e^ Pl ) hl . . . (e^p L ) h -) = s A (y). 

The calculation for G_£,/ 2 (t) is similar. Using the Cauchy-Littlewood identity ( 13 .3(1 . we 
get: 



G±L/ 2 (t) 



w({Mi» Mi!... 



exp f T 7 E^ + EE fcrfe^) . (3.16) 

V n=l fc>la=l / 



Let \u) = \u({Mi})) e U({Mi}) be an eigenstate of T(u,t), 

T(u,t) \u({Mi})) = \u({Mi})) r u (t; |u;({M t })», 
then the corresponding eigenvalue r u (t; |o;({Mj}))) can be written in the form 

L 

r„(t; |w({Mj}))) = C(t) J] sinh( 7 ( M - Ufc (t)). (3.17) 

fc=i 

We will call the expression in the right hand side a trigonometric polynomial (of degree 
L). The common multiplier C(t) and the roots of this trigonometric polynomial depend 
on all the times ti,t 2 , ■ ■ ■ (and on \oj)). Comparing with A3. 16j) . we find 

2 L D N 

c{t) = m 1 \...m n \ ^(EDtfi^W))- 

From (13.71) and (12.111) it is clear that the initial values of these roots are inhomogeneity 
parameters at the lattice sites: ttj(O) = U{. 



4 Trigonometric solutions of the MKP hierarchy 

In this section we study solutions of the MKP hierarchy which are periodic in the variable 
t = u with period 2717/7. We call them trigonometric solutions. For solutions of this 
class, the r-function is a "trigonometric quasi-polynomial" of u, i.e., a Laurent polynomial 
of the variable e 7 " possibly multiplied by an exponential function of u. 



4.1 The construction of trigonometric solutions 

By trigonometric solutions of the MKP hierarchy we mean r-functions which are polyno- 
mials in e 7 " for some 7 multiplied by an exponential function of u. They can be viewed as 
degenerations of double-periodic (elliptic) solutions in the complex plane of the variable 
u = t (they correspond to vertex models with elliptic R- matrices) . The general theory 
of elliptic solutions for the KP hierarchy was developed in [22] and extended to the MKP 
hierarchy in [26J. The trigonometric degeneration simplifies the construction and makes 
it more explicit [27]. Here we apply it to the case of our interest. The r-function for the 
trigonometric solutions will be obtained below in the form of the Casorati determinant 



Let r u (t) be the r-function of the MKP hierarchy. The Baker- Akhiezer function and 
its adjoint are defined in the following way [19]: 



^(M) = z"e^> T " (t .|f 1]) , (4.1) 

r, t 



<(t, z) = z~ u e-^ Tu[t+ } Z 1]) ( ! .2 ) 



T u (t 
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(£(t, z) is given in (j3.10j) ). In general, the ratios of the r-functions in the right hand sides 
can be expanded in infinite series around z = oo: 



WM) = rf,. 1 + + + , (4 . 3) 



«( M ) = 2-»e-«Mpi^ + ^ + ...V (4.4) 

According to the Krichever's theory of general algebro-geometric solutions [21], these 
solutions can be characterized and explicitly constructed by fixing certain analytic prop- 
erties of the Baker-Akhiezer function on a Riemann surface of the complex variable z 
(the classical spectral parameter). Recall that the quantum spectral parameter u is the 
zero time in the classical MKP hierarchy. 

For trigonometric solutions, the Riemann surface is the Riemann sphere (compact- 
ified complex plane), which represents a genus zero algebraic curve with singularities. 
Correspondingly, the Baker-Akhiezer function is, in this case, a rational function on the 
complex z-plane multiplied by power-like and exponential factors which give the required 
asymptotics (the essential singularity at infinity). For non- integer u the points and oo 
are branch points for the the Baker-Akhiezer function. In order to make it single- valued, 
one should make a cut between and oo. 

We know that the second series in (13. lip truncates at a = N. This suggests to 
assume the following ansatz for the Baker-Akhiezer function, in which the series in (14. 3 p 
truncates at the iV-th term: 

WM)=z VM(l + ^ + ... + !^!). (4.5) 

This explicitly defines the function z^e'^^'^ip^t, z) as a rational function on the ex- 
tended complex plane. The multiplicity N of the pole at z = is a discrete parameter 
characterizing the class of solutions to be constructed. Fix iV points pi G C, N non- 
negative integer numbers Mj such that 



M x + M 2 + . . . + M 



N 



and the set of parameters b i rn with i = 1, . . . , N, m = — |Mj, |Mj + 1, . . . , \M, t 
assume that b i ±i M . ^ for all i). Let us impose iV conditions of the form 



we 



£ b i>n ^ v (t,p i e 3 ^ m ) = 0, i = l,...,N, (4.6) 

rn=-Mi/2 

which are supposed to hold for any values of u, t». The sum goes over all integer numbers 
between — \Mi and —\Mi for even Mi and over all half- integer numbers between — \M{ 
and -\Mi for odd M». 

These conditions yield a system of iV linear equations for iV coefficients Wk which 
allows one to fix the Baker-Akhiezer function ip. The general theory guaranties that the r- 
function associated with this -^-function according to ( 14.ip solves the MKP hierarchy. The 
points pi and entries of the matrix b^ m are parameters of the solution. The coefficients Wk 
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appear to be rational functions of e 7U while the r-function is a trigonometric polynomial 
in u (possibly multiplied by an exponential function of u). From the algebro-geometric 
point of view, these solutions are associated with singular Riemann surfaces with N 
"strings" of singular points 

with the center at p;. The points of each string are glued with each other in a complicated 
way. Note that the parameters bi >m can be multiplied by any non-zero complex numbers 
kf. the transformation bi )Tn — > kibi >m does not change anything. 

The family of periodic iV-soliton solutions is a very particular case of this con- 
struction corresponding to Mj = 1 for all i. In this case conditions (14. 6p become 
bi,-i/2^u (t, Pie -7 ) = —bi^^ipui^^i^ 1 ) and the solutions are associated with the Riemann 
sphere with N pairs of double points p^e 7 and p^e -7 . 

It is easy to see that conditions ( 14 .6p are equivalent to the system of linear equations 



N 



Bi{u, t) + B i(u - k, t)w k = 0, 



(4.7) 



k=l 



where 



bi, m {t) 



+Mi/2 

t):=p? J2 b ijm (t)e 2 ^ mu 

m=-Mi/2 

The system can be solved using the Cramer's rule. This gives the following explicit 
expression for the Baker- Akhiezer function: 



1 z- 1 
Bi(«,t) B x (tt-l,t) 



-A' 



Bt(u-N,t) 



B N (u,t) B N (u-l,t) ... B N (u-N,t) 



B x (u -l,t) 



BUu-N^) 



B N (u -l,t) ... B N (u-N,t) 
Comparing with (14. ip . we conclude, using the obvious property 

Bi{u, t - [z- 1 ]) = Bi{u, t) -Bi(u+ 1, t)*" 1 , 



(4.9) 



(4.10) 



that the r-function is given by the difference Wronskian (Casorati) determinant in the 
denominator: 

r«(t) = (detg)- u . det B^u-j, t) (4.11) 

i,j=l,.. .,N 

(here g = diag(pi, . . . ,pjv) is the same matrix as in the previous sections). The factor 
(detg) - " is put here to make r u (t) a pure trigonometric polynomial in u of the form 
(13.17p . Comparing the highest and the lowest coefficients in (I3.15P (given by equation 
(I3.16P ) with the corresponding coefficients in (14. lip , we see that the parameters b i± \ M _ 
obey the following relations: 



N 



i=l 



±lM k 



j<k 



(±1) L L\ 



M 1 \...M J 



EN 
n=l "» 



N ■ 
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From ( 14. 9 p it is clear that the last coefficient in (14. 5p . w^, in terms of the r-function 
is given by 

w N (u,t) = 



V T u+l (t) 



T u (t) • 



We also note the formula 



u;i(w > t) = -3 Cl logT tt (t) (4.12) 
for the first coefficient in (14.51) . w±, which easily follows from the obvious relation 

d tl B i (u,t) = B i (u + l,t). (4.13) 



Rewriting (I4.10p in the form 

Bi(u, t + [z- 1 ]) = Bi(u, t) + z- x Bi(u + 1, t + [z- 1 ]), 
it is straightforward to check that 



r u (t + [z- 1 ]) = (detg)- u 



B 1 {u-l,t+[z- 1 }) Bi(«-2,t) 
B 2 (u-l,t+[z- 1 ]) B 2 {u-2,t) 



B N {u-l,t + [z- 1 ]) B N (u-2,t) 
It directly follows from the definition that 



Ei(w-iV,t) 
B 2 (u-N,t) 

B N {u-N,t). 



(4.14) 



M i/ 2 i, 

BAu, t + [Z- 1 ]) = V u t V 0i ' m * e ^muH(t, Pi e^)_ 

m=—Mi/2 

Expanding this in powers of we get: 

Bi(u, t + [z- 1 ]) = Bi(u, t) + Bi{u + 1, t)z _1 + Bi(u + 2, t)z~ 2 + 
Therefore, the expansion of r u (t + [z^ 1 ]) around oo reads 



(4.15) 



l (t + [z- 1 ]) = (det^)-^ 



s=0 



B l (u+s-l,t) Si(u-2,t) ... Bx(u-N,t) 
B 2 (u+s-l,t) B 2 (u-2,t) ... A 2 (u—N, t) 

B N (u+s-l,t) B N (u-2,t) ... B N (u-N,t) 



(4-16) 

We thus see that the adjoint Baker- Akhiezer function has the determinant representation 

B x (u-l^ + lz- 1 ]) £i(u-2,t) ... Si(u-JV,t) 
5 2 ( M -l,t + [^ 1 ]) B 2 (u-2,t) ... S 2 (u-JV,t) 

J B 7V (M-l,t + [^ 1 ]) Bjv(«-2,t) ... B N (u-N,t) 



r u (t,z) = z- u e 



Si(«-l,t) £i(u-2,t) ... 5i(u-iV,t) 
S 2 («-l,t) S 2 («-2,t) ... B 2 (u-iV,t) 

Sjv(u-l,t) S w (u-2,t) ... B N (u-N,t) 



(4.17) 
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Let us introduce the notation 



B k (u, t) := det Bi{u + 1 - j, t) (4.18) 

i=l, ...,^,...,iV 
2=1, ... , N-l 

for the minor M k ^ of the NxN matrix Bi(u + 1 — j), 1 < z, j < N. Then, expanding 
the determinant in the numerator of (14. 17ft in the first column, we obtain: 

-u -£(t,z) n 

C(t, z) = — J2(-l) k - l B k (u - 2, t) B k (u - 1, t + [z- 1 ]), 

or, substituting (I4.15p . 

(det^r u (t)^ 1 m jf-. /a ^- Pl e 2 ^ 

(4-19) 

This gives the pole expansion of the adjoint Baker-Akhiezer function. We see that in 
general it has simple poles at all the points forming the "strings". Below we need this 
formula rewritten for the function r u (t + [z -1 ]): 



N Mi/2 , 27m(u-l)+^(t,pie 2 T m ) 

r u (t + [z- 1 ]) = z(detg)- u Y,(-n- 1 Pr 1 E — 2^ A(«-2,t), 

i=l m=—Mi/2 

(4.20) 

with simple poles at the same points. 



4.2 Undressing Backlund transformations for the trigonometric 
solutions 

As it was demonstrated in [1] for models with rational -R-matrices, the main relations 
of the Bethe ansatz method are naturally built in the construction of rational solutions 
to the MKP hierarchy. The nested Bethe ansatz scheme appears to be equivalent to a 
chain of some special Backlund transformations of the initial rational MKP solution that 
"undress" it to the trivial solution by reducing the number of singular points in succession. 
All this remains valid for vertex models with trigonometric .R-matrices, with the only 
difference that the undressing procedure should be applied to the trigonometric solutions. 
Technically it becomes even simpler because poles of the adjoint Baker-Akhiezer function 
are simple in this case. In particular, the functions B k (u, t = 0) and B k (u, t = 0) should 
be identified, up to some irrelevant factors, with the (eigenvalues of) the Baxter Q- 
operators on the first and the last levels of nesting in the nested Bethe ansatz scheme. 

Adding or removing a "string" pie 2im with the center at pi to or from the data of a 
trigonometric solution is a Backlund transformation. It sends a trigonometric r-function 
to another one. We will be interested in the removing of a string that results in decreasing 
the degree of the trigonometric polynomial (the undressing transformations). Basically, 
such a transformation can be done by extracting the singular part (the residue) of the 
function r u (t + [z^ 1 ]) at any of its simple poles which are located at the points of the 
string with the center at Pi. Specifically, consider the function 

r u W(t) = (-l) i - 1 pr 1 det^e^( tl+1 )-«( t ^ e ^" I )res 2=Pie - 7 M I r u+1 (t+ [z^ 1 ]). (4.21) 
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Equation (I4.20p implies that 

T^\t) = b^ ¥k {p- 1 detg)- u B l {u-l ) t) 1 

i.e., up to the irrelevant constant factor, it has exactly the same determinant form as 
T u (t) with the string with the center at pi removed. Therefore, it is a r-function, i.e., 
it satisfies the same Hirota equations as r u (t) does and r — > is indeed a Backhand 
transformation. The degree of the trigonometric polynomial rj*'(t) is L — Mj. Note that 
the residue in f 14 . 2 lj) is taken at the left edge of the string. This has an advantage that 
the coefficient b i _i M . is non-zero by definition and thus the result of the transformation 
never vanishes identically (the same holds for the right edge). 

The procedure can be continued until one obtains a polynomial of degree 0. The 
inductive definition is as follows. Fix a set I n = {ii, . . . , i n } C {1,2,..., N}. Suppose 
we have a r-function r^ ll2 "' ln -^(t) obtained at the (n — l)-th step, then the r-function 
at the n-th step is defined as 



(4.22) 



r [iifa...Hd(t) = (-1)--^ n Pi)(^,-j«j" i e 7W< - (,H - i) - €(t ^ e " 7 in) 

je{i,...,iV}\j„ 

xres^^r^-^t+K 1 ]). 
This function has the determinant representation 

r M,..*n] (t) = ^ 6 ._ |m J( n pJ «) ^ Bi{u _ j} t) . {m 

i€ln j£{l,...,N}\I n j=l,...,N— n 

As it is shown in detail in pQ , the "undressing" chain of Backlund transformations 

r u (t) r^(t) r^^(t) . . . r|-^(t) 

is equivalent to the nested Bethe ansatz scheme, with the r-functions l *J (t) being 
eigenvalues of the master T-operators on higher levels of the nesting procedure. In 
particular, r| n """ l "l(t) at t = are eigenvalues of the Baxter's Q-operators. They are 
trigonometric polynomials in u of decreasing degree as n increases. This implies the 
system of Bethe equations for their zeros. 



5 Zeros of the master T-operator as the Ruijsenaars- 
Schneider particles 

As we have seen, eigenvalues of the master T-operator are trigonometric polynomials in 
the spectral parameter u of the form (I3.17p . The roots of each eigenvalue have their own 
dynamics in the times £j. This dynamics is known [26] to be given by the trigonometric 
Ruijsenaars-Schneider model [6]. The inhomogeneity parameters Ui are coordinates of 
the Ruijsenaars-Schneider particles at U = 0: Ui = Ui(0). 

Here we derive, following [26], the equations of motion for zeros of the trigonometric 
r-function (the master T-operator) with respect to the first time flow t% = t. Our starting 
point is the differential-difference equation for the Baker- Akhiezer function: 

d tl ip u (t, z) = Vv+-i(t, z) + d tl \og Tu+1 /^ ip u (t, z). (5.1) 
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which follows from the definition and from the Hirota equations. 

It is clear from (14. ip that ijj u {t) has simple poles at u = Uj(t). Let us introduce the 
function 

sinh(7(-u + 0) , x -, / >n 

$ M > = ■ r / ; ■ w = coth( 7 n) + coth( 7 C) 
smh(7ttj smh(7Q 

and adopt the following pole ansatz for the Baker- Akhiezer function: 

L 

^(t, z)=z u J2 Si(t, 2, C))$(« - %(t, 0- (5-2) 

Here £ plays the role of an auxiliary spectral parameter. Substituting this ansatz into 
(15. ip . one is able to derive the equations of motion together with their Lax representation. 
Skipping further details of the calculations, we give the results. The double poles at 
u = Uj cancel automatically. Cancelation of simple poles at u — Uj — 1 yields: 

1 u j Yl ~ X > s k = zsj, j = 1, . . . , L, 

k 

where Uj := d tl Uj, Uj k := Uj — u k . Cancelation of simple poles at u = Uj yields: 
= l u j J2 ®( u jk, s k + 7[coth(7C) Uj + J2 u k (cothj(Uj k ) - coth^u^+l)) 



Finally, comparison of the constant terms at u — > ±oo (we assume that 7 is real positive) 
yields the condition 

X) h = z X s 3 

which does not add any new constraint because in fact follows from the previously ob- 
tained relations. The conditions obtained above can be written in the matrix form as 

£(C)s = zs, s = M(C)s, (5.3) 

where s = (si, s 2 , . . . , sl)* an d the matrices C((), Ai(C) are defined as 

^(0 = 7^^(^-1,0, (5-4) 

M?fc(0 = 7 (coth(7C)-coth7)w j + ^ V(uji) 8 jk + 7(1 - 5j k )uj$(uj k , (), (5.5) 

where 

V(u) := coth( 7 w) — coth(7(u + 1)). 
The compatibility of equations (I5.3P implies the Lax equation 

£(C) = [M(0, C(()}. (5.6) 

A direct calculation shows that it is equivalent to the equations of motion for the Ruij- 
senaars-Schneider system: 



u 



.1 



^Uj^u k (y(u jk )-V(u k , 



k+j 



(5.7) 

-2 7 sinh 2 7 Y ^ fc cosh( 7%fc ) 

^ sinh(7(u jfc -l)) sinh(7M jfc ) sinh(7(w ifc + l)) 
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In the course of the calculation, the following identities are useful: 



$(w - 1, c)$(v, C) - 0$(« - 1, = $(« + v - l, C) - y(- 



-1,0 = 7(coth( 7 C) - coth 7 - j$(ti -1,0- 7$(-l, 0$K 0- 

We also note that the system (15.61) is a Hamiltonian system with the Hamiltonian 

Ur = f & n ( sin M7(^ + l))sinh(7K. fc -l))\ 1/2 
ft £y\ sinh( 7 u jfc )) ; 

and the canonically conjugate variables with the Poisson brackets {vj,Uk} = 5jk- 

There are also higher Hamiltonians Uj in involution which generate the higher flows with 
respect to tj. 

The spectral curve is given by the equation 

det (£(C) - z) = 0. (5.9) 

One can show that this curve is the Riemann sphere with points of each string pie 2 " 1 ™ 
being glued in a complicated way. The coefficients of the characteristic polynomial in 
the l.h.s. are integrals of motion for the Ruijsenaars-Schneider system. 

Finally, let us stress the specific way of posing the problem in the context of the 
Ruijsenaars-Schneider system that corresponds to solution of the vertex model or quan- 
tum spin chain. The standard mechanical problem is: given initial coordinates and 
velocities of the particles %(0), %(0), find the time evolution Uj(t). By contrast, in order 
to find eigenvalues of the transfer matrix, one should pose the problem in the following 
non-standard way: given initial coordinates Uj = Uj(0) and values of all higher integrals 
of motion Uj, find initial velocities Uj(0)- Indeed, the initial velocities allow one to restore 
the transfer matrix via residues at its poles: 



T{u) 



T%u,0) 



= -Mo). (5.10) 

H{Mi})) 



The solution is not unique: different possible solutions to this problem correspond to 
different eigenstates of the transfer matrix in the sector "H({Mj}). 
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